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Abstract 

A simple formalism describing the light response of CsI(Tl) to heavy ions, which 
Preprint submitted to Elsevier Preprint 7 February 2008 



quantifies the luminescence and the quenching in terms of the competition between 
radiative transitions following the carrier trapping at the Tl activator sites and the 
electron-hole recombination, is proposed. The effect of the 5 rays on the scintillation 
efficiency is for the first time quantitatively included in a fully consistent way. The 
light output expression depends on four parameters determined by a procedure of 
global fit to experimental data. 

Key words: PACS number: 29.40.Mc, 32.50.+d 

(light response of CsI(Tl) to heavy ions, quenching, delta rays) 



1 Introduction 

The thallium-activated caesium iodide (CsI(Tl)) scintillation crystals combine 
the advantages of high stopping power and high reliability for large thickness, 
good energy resolution (a few percent) and the possibihty of hght particle 
isotopic identification using the pulse shape discrimination technique. The 
response of these detectors has a non-linear dependence on the energy of the 
incident particle. The reduction of the light yield with regard to the linear 
behaviour is referred to as "quenching". In addition, for a given energy, the 
light output depends on the type of the particle. 

Obtaining a particle dependent light output description and energy calibra- 
tion of these scintillators, over a very large range of energies and ion atomic 
numbers, is the primary objective of the present work. Quite successful early 
approaches of the main experimental observations to be considered are dis- 
cussed in section 2. A light response expression depending on the mass, charge 
and energy of the ionizing particle is derived in section 3, in a simple formal- 
ism which makes use of band theory scheme in insulators. The quenching is 
connected to the carrier recombination or their trapping at crystal imper- 
fection sites, eventually enhanced by the nuclear interaction of the slowing 
down fragment with lattice nuclei. The contribution of the generated 5 - rays 
is quantitatively included. The capability of this "recombination and nuclear 
quenching model" (RNQM) to describe experimental data is demonstrated in 
section 4. These data were obtained by means of the INDRA array, through 
a procedure described in the accompanying paper [1]. A summary is given in 
section 5. 
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Notation and values of physical constants and variables used in this paper. 



Symbol* Definition Units or Value 

Physical constants 

c speed of light in vacuum 299 792 458 m s"^ 

u unified atomic mass unit 931.494 32 MeV/c^ 

rUe electron mass 0.510 999 06 MeV/c^ 

AfA Avogadro constant 6.022 136 7x 10^^ mol"^ 

k Boltzmann constant 8.617 385 x 10"^ eV K"^ 

h Planck constant 41.356 692 x lO'^^ j^eV s 

CsI(Tl) crystal 

p Csl crystal density 4.53 g cm~^ 

< Acsi > average atomic mass number for Csl 129.905 

< Zcsi > average atomic number 54 

Nh host cation and anion volume ions cm~^ 
concentrations 

Nao Thallium activator volume concentration ions cm~^ 

C Thallium molar concentration % 

e energy gap between the valence and con- 10 eV [2] 

duction bands 

Charged particle passing through the CsI(Tl) crystal 
A atomic mass number 
Z atomic number 

E current kinetic energy in the CsI(Tl) MeV 
crystal 

X coordinate along the path 10^^ atoms cm~^ 

Eq experimental initial kinetic energy in the MeV 
crystal 

R{Eq) range 10^^ atoms cm~^ 

Qo approximate total integrated charge oc ex- a.u. 
perimental light output 

a, (3 relativistic kinematic variables 

eg energy per nucleon threshold for 5 - ray MeV/n 
production 

Eg energy threshold for S - ray production MeV 
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Symbol* 


Definition 


Units or Value 






the corresponding relativistic kinematic 
variables 








AE 


energy deposited in the preceding detec- 
tion layer 

5 - rays 


MeV 






T 


kinetic energy transferred to an electron 








T 

max 


maximum kinetic energy transferable to 
an electron 








Tcwt 


value of T above which an electron is con- 
sidered as a (5 - ray 








Re 


electron range 


mg cm~^ 








radius of the highly ionized primary 
column 

Stopping powers 


cm 






A 


density effect in Bethe-Bloch formula 








7 . . 


effective charge of the impinging particle 








7e// 


coemcient oi enective cnarge 








r 
J 


ionization constant lor KjSL 


579.8 eV 






— ur!j/ ax 


specific energy loss 


eVx 10-^5 


atoms" 


-1 2 

cm 


-(dE/dx)e = Se 


specific electronic stopping power 


eVxlO-^^ 


atoms" 


-1 2 

^ cm^ 


-(dE/dx)n = Sn 


specific nuclear stopping power 


eVxlO-i^ 


atoms" 


-1 2 

cm 




electronic or nuclear stopping power 


McV cm"^ 




-(dE/dx)T<Tcut 


restricted specific electronic stopping 
power 


eVx 10-15 


atoms" 


-1 2 

cm 


-(dE/dx)p 


part of the specific electronic energy loss 
remaining inside the primary column 


eVx 10-15 


atoms" 


1 9 

cm 


-(dE/dx)5 


part of the specific electronic energy loss 
carried outside the primary column by the 

O 1 ciy o 


eVx 10-15 


atoms" 


^ cm^ 




fractional energy loss transferred to a 5 - 










TTi Pfi n pn prp"v ypn 1 1 1 rpH to nrorl 1 1 rp a 
dislocation 


eV 






Nn 


concentration of anion or cation defects 
produced by nuclear interaction 


ions cm-^ 






e' 


cm. deflection angle 


rad 






da'/ an' 


cm. Rutherford cross section 


cm 






^Ruth 


concentration of aB^on or cation defects 
produced by Coulomb scattering 


ions cm"^ 







Symbol* 


Definition 


Units or Value 




Calculated light output and related variables 


L 


Ccilculcited integral light output 


a.u. 


S 


oCillLlilclLlOIl cliiClcilCy idLLUl Ui IDnKb 


a.u. 


JCB 


quenching factor of Birks 


MeV"^ cm^ g-^ 




infinitesimal light output inside the pri- 
mary column 


a.u. 


(dL/dx)^ 


infinitesimal light output outside the pri- 
mary column 


a.u. 


hu 


bCiiiLinaijion eneigy 


e V 


Wd 


energy of a nearest-neighbour bound in 
solid 


« 1 eV 


Nd 


concentration of anion or cation defects 
produced by thermal vibrations 


ions cm"^ 


T 


absolute temperature 


°K 


t 


time 


s 


rie 


current volume concentration of free elec- 
trons created by the ionizing particle 


electrons cm~^ 


nh 


current volume concentration of holes cre- 
ated by the ionizing particle 


holes cm~^ 


no = neo = Uho 


initial volume concentration of carriers 
created by the ionizing particle 


electrons cm~^ 


Na 


current volume concentration of electron 
traps at Tl activator sites 


ions cm~^ 


ANa 


current volume concentration of hole traps 
at Tl activator sites 


ions cm^'^ 


AAe 


coefficient of electron trapping probability 
at activator sites 


s~^ cm^ 


^Ah 


coefficient of hole trapping probability at 
activator sites 


s cm'^ 


^De 


coefficient of electron trapping probability 

at defect sites 


s~^ cm^ 


^Dh 


coefficient of hole trapping probability at 
defect sites 


— 1 'H 

s cm 




coefficient of trapping probability at de- 
feet sites for both type of carriers 


s cm'^ 




coefficient of carrier recombination 
probability 


s"^ cm^ 


A' 


AAe,D,R/(AAe^A + AdATd) 


cm^ 


■^calc 


calculated initial kinetic energy in the 
CsI(Tl) crystal ^ 


MeV 


a 


multiplicative constant 


a.u. 



Symbol* 


Definition 


Units or Value 




gain fit parameter in the exact expression 


a.u. 




nf T 




aR 


recombination quenching fit parameter in 


eV~^10^^atoms~^cm^ 




the exact expression of L 




an 


nuclear quenching fit parameter in the ex- 


eV~^ 10^^atoms~^cm^ 




act expression of L 





*Most of the notations of the original references have been kept. 



2 The total light response of CsI(Tl): two historical hints 

2. 1 Quenching 

The CsI(Tl) light response to an ionizing charged particle of atomic number 
Z and mass number A is non-linear with respect to the incident energy Eq, 
depending on Z and A [3-5]. The quenching has been associated with the 
specific energy loss of the particle. This property provides a means for the 
pulse shape discrimination technique. The differential light output per unit 
path length expression dL/ dx = S\ dE/ dx|c/ (1 + 1CB\ dE/ dx|c), proposed by 
Birks for alpha particles in organic scintillators [6] and applied also, without 
rigorous support, to inorganic ones [7], leads to the differential scintillation 
efficiency: 

dE l + /C-B|dE/dx|e' ^ ' 

where S and ICB are the scintillation efficiency and the quenching factors, 
respectively [7]. The differential scintillation efficiency is monotonically de- 
creasing with specific electronic stopping power — ( dE/ dx)c. This relation was 
applied with reasonable results in case of light charged particles or intermedi- 
ate mass fragments {Z < 15) [7,8], i.e. as long as the fraction of the light yield 
from (5-rays is not significant. Under the approximation | dE/ dx|e oc AZ^/E, 
Eq. (1) was integrated and successfully used to describe the total light output 
from fight ions [9]. 

2.2 Knock-on electrons 

Experimental studies, performed with ions of a few MeV/nucleon, as heavy as 
those from Boron to Neon [10,11], have shown that the scintillation efficiency 
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Fig. 1. Experimental light efficiency versus the estimate of the average specific elec- 
tronic stopping power, for different fragments (AZ = 5). Data from Xe + Sn system 
at 32 and 50 AMeV [1]. The discrete loci are the effect of the 5 - rays. Its description 
in the framework of the present scintillation model - solid lines - is compared to the 
experimental data in the upper box (AZ = 10). 

of alkali halide crystals (Nal(Tl), CsI(Tl)) is not a function of — (dE/ dx)e 
alone, but instead, is composed of a series of discrete functions, one for each 
incident particle. The observation has been thoroughly analyzed by Murray 
and Meyer [12], as shown in this subsection. 

This effect is illustrated with INDRA data from Ref. [1] in Fig. 1, where the 
experimental light efficiency {Qq/Eq) is shown, for different heavy ions, versus 
the estimate of the average specific electronic stopping power AZ'^/Eq [1]. Qq 
is the experimental light response corresponding to Eq. There is a domain 
where, at the same average specific energy loss, the higher the atomic num- 
ber Z, the higher the scintillation efficiency. The above mentioned discrete 
functions have been interpreted as a signature of modestly energetic (~ 1 
keV) knock-on electrons, the alternative term for 6 - rays, whose ranges ex- 
ceed the radius of the high carrier concentration column along the wake 
of the ionizing particle, estimated to be around 40 nm [12]. The differential 
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light output per unit path length is the sum of two contributions: one from 
the highly ionized primary column and one from the knock-on electrons. Con- 
sidering the total specific electronic stopping power of the incident particle: 
( — dE/ dx)e = ( — dE/ dx)p + ( — dE/ dx)^ as the sum of the energy deposited 
inside the primary column (— dE/ dx)p and the energy carried off by the S 
- rays (— dE/ dx)^, the fractional energy loss JF of the incident particle, de- 
posited outside the primary column, was defined [12]: 

(-dE/dx)^ 

(-dE/dx)e- ^ ' 



The differential scintillation efficiency is then given by: 



The scintillation efficiency referring to the primary column (dL/dE)p may 
be treated with formula (1) provided that (— dE/dx)p is known. The 5 - 
ray scintillation efficiency ( dL/ dE)^ is very nearly a constant for energies 
higher than 1 keV [2,13]. Estimations of T have shown that this quantity 
depends only on the energy per nucleon E/Aoi the particle and it is different 
from zero only above a certain threshold value eg — (-B/^) threshold [12]- These 
characteristics will be found again in the next section in a very simple approach 
for calculation of valid for all ions. 

Thus, for a particle slowing down in a CsI(Tl) detector from the incident 
velocity oc Eq/A > ^Jel to zero, the differential scintillation efficiency is 
given by Eq. (3) as long as (E/A) > e^, and by Eq. (1) in the last part of 
its range, when (E/A) < es- The total light output would be obtained by 
integrating the scintillation efficiency over the entire range of the particle. 
One has to note that the above formulae (1), (3) are not time dependent and 
consequently they cannot predict the time behaviour of a scintillator signal. In 
the next section, we shall develop a simple time dependent formalism leading 
to scintillation efficiency expressions which, up to a supplementary nuclear 
quenching term, reduces in the first order approximation, to relations very 
similar to equations (1), (3), when integrated over time. 



3 Recombination and nucleeir quenching model 
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3.1 Experimental evidence 



The Csl scintillator is an alkali halide crystal having a body centered cubic 
structure. The present work concerns data from CsI(Tl) operated at room 
temperature and having a given molar concentration C ~ 0.02% — 0.2% [14]. 
The corresponding thallium volume concentration is A^^o ~ 2.11 x 10^^ — 
2.10 X 10^^ ions cm~^. Experimental observations related to the scintillation 
and its behaviour as a function of temperature in both pure Csl and activated 
CsI(Tl) crystals can be found in Refs. [7,15-25]. In the latter case, the thaUium 
concentration dependence is regarded too. At room temperature, the emission 
spectrum of a Csl(Tl) excited by charged projectiles [16] is dominated by a 
broad yellow band, centered near 550 nm. This band is thus associated with 
emission from Tl luminescence centres and is referred to as the Tl band [16]. 
The yellow band nearly completely replaces the competing blue band, related 
to /~ vacancies or elementary colour centres [22,23] and the ultraviolet (UV) 
band, considered as an intrinsic property of a perfect crystal [16]. Both of the 
latter bands are present in case of pure Csl. 

We briefly remind the experimental evidence related to CsI(Tl) at a given 
temperature and concentration, on which our model is dependent, a) The 
scintillation efficiency is non-linearly increasing with the increasing energy of 
the particle. For different ions at the same energy, the scintillation efficiency 
decreases with increasing specific electronic stopping power [7] . If the ionizing 
particle is an electron or a 7 - ray of energy higher than 100 keV, the light 
output may be considered to a good approximation as a linear function of 
energy [24,26]. b) The scintillation efficiency depends on the fraction of de- 
posited energy creating 5 - rays [12], which affects the monotonic behaviour 
resulting from the previous item, c) Combined studies of light efficiency and 
scintillator decay-constant versus temperature [7,27] have indicated that the 
quenching effect is external to the luminescence centre [7]. 



3.2 The scenario 

An ionizing particle excites the scintillator crystal along its path by creating 
carriers: electrons, promoted in the conduction band, and associated holes in 
the valence band. CsI(Tl) is an insulator and the energy gap between these 
two bands is of the order of the mean ionization energy e 10 eV [2] . ThaUium 
iodide is melted with the caesium iodide host bulk before the crystal growth 
[14]. Some cation sites are occupied by Tl"*" ions having lost one valence electron 
in favour of I^ anions. As in the case of KCl(Tl) [28], the excited energy levels 
of the T1+ ions are presumably localized in the forbidden band of the host 
lattice. Even though deeper than in the semiconductor case, these levels may 



9 



essentially participate in the ionized lattice deexcitation, provided that they 
are free. 

Although the emission bands of the CsI(Tl) scintillators, of which the most 
prevalent one is centred at about 550 nm, have been intensively studied, 
there is still controversial discussion about the origin of this band [28-40]. 
Thus, according to different interpretations, the luminescence emission may 
imply, for example, a single thallium centre [31], a couple of T1+ - Vk centres 
[38] or even more complex cluster configurations nearby a Tl impurity [39,40]. 
As definitions of the imperfection centres we use those in Ref. [41]. Most of 
the mentioned approaches converge towards two ideas: for the luminescence 
induced by highly ionizing particles, both type of carriers, electrons and holes, 
play a role in the emitted light, i.e. in the scintillation efficiency; the radiative 
emission around 550 nm takes place in the {Tl~^)* excited ions. Let us imagine 
a formal, oversimplified scenario for the scintillation process, by exemplifying 
it on the latter luminescence centres. 

We suppose that the T1+ ion, with its 2 remaining valence electrons, super- 
fluous in the ionic bond of the Csl lattice, has a donor behaviour. Stimulated 
by thermal vibrations, the Tl"*" ions are loosing, in part, some of the outer 
electrons, creating intruded levels into the upper part of the forbidden band. 
These double ionized thallium atoms Tl"*""*" would become electron traps of 
volume concentration Na- The remaining T1+ ions would be the hole traps, 
with the volume concentration AA^^ = Nao — Na- One could reason as fol- 
lowing: in the highly ionized fiducial volume, the holes are mainly trapped, 
via a mechanism which is disregarded here, by T1+ ions becoming Tl"*""*" ions, 
while the electrons are mainly trapped on excited levels by such newly created 
double ionized thallium atoms or by previously existing ones due to thermal 
excitation. They form (TZ^)* excited ions which decay afterwards by emitting 
the scintillation of energy hu, in the characteristic 550 nm band; here h is the 
Planck constant and u is the light frequency. The scheme of the process is: 

ri++ + electron — ^ (T/+)* + hu 
Tl+ + hole — >Tl++. 

The mathematical formalism is the same if a more complex configuration of 
the luminescence centre is considered. During the non-equilibrium state in the 
cylinder around the ionizing particle trajectory, the two states of thallium feed 
one another via the carrier trapping, closing in this way the cycle which insures 
the restoration of the equilibrium and the local electrical neutrality. A time 
dependent formalism which treats the kinetics of the volume concentrations 
for both types of carriers and thallium traps would somehow cover also the 
"exciton" scenario [2,7,37-40]. In this way, it would not be necessary to isolate 
in time the two steps of the scintillation process: prompt formation of exciton 
and its delayed deexcitation. 
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Similarly to the levels created by thallium impurities, deeper levels may be 
intruded by other point defects of the crystal: the cations and anions removed 
towards interstitial positions, or for the pure alkali halides more probably to 
the crystal surface [21] , and the associated vacancies. The positive and negative 
defects are produced by thermal vibrations in practically equal concentrations 
Nd [21]. These levels may compete with the activator sites in carrier trapping, 
following radiative transitions in the blue band or non-radiative transitions. In 
insulators, the Fermi-Dirac statistics of the carrier states is well approximated 
by a Boltzmann function. The concentration of those carriers located at the 
defect sites, together with the neutrahty condition [42], provide the defect 
distribution, described by a Boltzmann function too. The energy of a nearest- 
neighbour bound in a solid is of the order of 1 eV. Actually, after Kittel [21], the 
cation and anion defect concentrations are Nd pa x exp{—WD/kT), where 
Nff is the host volume concentration of cations and anions and Wd ~ leV 
is much larger than the mean vibrational energy. Comparable values for Wo 
(0.2 - 0.4 eV) are obtained with the above expression from concentration of 
elementary colour centres of the order of 10^^ - 10^^ centres cm^'^, determined 
by means of absorption coefficients [41]. In a Boltzmann type distribution, Njj 
would represent the number of cations or anions per unit volume having ener- 
gies higher than Wo- Consequently, Wo ~ 1 eV would be the minimum energy, 
transferred to a lattice ion, required to produce a defect. At a given tempera- 
ture, the number of defects may be enhanced by the nuclear interaction of the 
incident particle with the lattice nuclei. This interaction is proportional to the 
specific nuclear stopping power — (dE/ dx)n = Sn and produces an additional 
defect concentration (practically the same for cations and anions) N„. 

Besides, one has to consider the capture of electrons by the ionized atoms in 
the wake of the particles, i.e. the "direct" electron-hole recombination, leading 
to the crystal intrinsic transitions in the UV band. This process probability 
increases with the specific electronic stopping power — (dE/ dx)e = Sg of the 
incident particle, while it has to be practically null for electrons. In experi- 
ments, a small UV peak was observed for the heaviest ions of low energy. On 
the contrary, for very energetic protons punching out through a thin CsI(Tl) 
crystal, in such a way that the Bragg peak is not contributing to the scintil- 
lation, both quantities Sg and Sn are so weak that the shape of their signal is 
nearly similar to that of electrons generated by 7 - rays. We do not consider 
the possible feeding of the yellow band by scintillations in UV band, a two 
step process in CsI(Tl) put in evidence by UV laser irradiation [43]. 

By quantifying these assumptions, it will be possible to account for the main 
experimental observations mentioned by items a) - c) of the subsection 3.1. 
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3.3 Stopping powers 



An incident fragment slowing in a CsI(Tl) crystal loses its energy mostly by 
ionization, and in a much smaller fraction, by interacting with the nuclei of 
the lattice. The specific energy loss is — ( dE/ dx) = Se + Sn = Se(l + Sn/Se), 
the sum of the specific electronic and nuclear stopping powers, expressed in 
eVxlO"^^ atoms"^ cm^. Multiplied by the crystal density p, they lead to 
the corresponding stopping powers per unit length: p x S^,^, expressed in 
MeVcm"^. For the calculation of the differential hght output, connected to 
the differential deposited energy, we have used the stopping powers of Ziegler 
[44], derived from master curves of a particles, based on Bcthc-Bloch formula. 
For the specific electronic stopping power cit energies above 2.5 AMeV, 
an effective charge Ze// — Z x 'yeff{E, A, Z) has been used with the pa- 
rameterization of Hubert et al. for jeff{E,A, Z) [45]. Below 2.5 AMeV, the 
parameterization of Ziegler was chosen for ^eff{E, A, Z) rcnormalizcd at 2.5 
AMeV, but not lower than 200 AkeV, where '~feff{E, A, Z) was kept constant, 
at the value reached at 200 AkeV: jeffiE = 200 AkeV, A, Z), for each type of 
fragment. At the atomic number Z, the integer of A = 2.072Z + 2.32 x IQ-^Z"^ 
for Z < 50 or A = 2.045Z + 3.57 x IQ-^Z'^ for Z > 50 [46] is used as mass 
number [1]. The resulting Se curves vs the energy per nucleon E/A are plotted 
- with solid lines - in Fig. 2a), for several atomic numbers. The corresponding 
specific nuclear stopping powers, as calculated by Ziegler's approach [44], are 
shown in Fig. 2b) - sohd fines. 

In order to determine J'{E) in the next subsection, wc introduce, besides the 
specific electronic stopping power (— dE/dx)e, the restricted mean rate of 
specific electronic energy loss (— dE/ dx)T<Tcut [48], which is the mean rate of 
energy deposited for collisions excluding energy transfers T to electrons greater 
than a cut off value Tcut- In this case, we have preferred analytical expressions 
of these quantities instead of the above mentioned Ziegler recipe which is used 
later for calculating light output and which is not easily to handle. These are 
based on Bethe-Bloch expressions that are valuable for moderately relativistic 
charged particles {E/A > SOMeV/u for the Csl medium), other than electrons 
[48], with the parametrization I — 16Z^^ for the ionization "constant" [49]. By 
grafting on these expressions the effective charge defined above Z^jf = Z^^ff, 
the lower limit of their domain of applicability goes down towards 1-2 MeV/u, 
as shown by the dashed curves in Fig. 2a) for ( — dE/ dx)e. 

The elastic Coulomb scattering on the crystal lattice nuclei - the main part of 

the nuclear interaction and hence of the specific nuclear stopping power - is 
used to quantify the density of lattice defects induced by the incident particle: 
A^Ruth = {p^fA/i{Acsl) X7rr2)) /; {da' /dn')dn' . Here (Acsi) = 129.905 g is the 

mm 

mean atom-gram, AOi is Avogadro's number and {da'/dfl') is the Rutherford 
cross section in the centre of mass (cm.) frame; the minimum deflection angle 
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in the cm. frame 



necessary to transfer an energy SE = Wd in one 
Coulomb deflection event, may be calculated in classical kinematics from the 
latter equality if the energy required to produce a defect in the lattice Wd is 
known. Then: 




E 



D 



{Acsi) 
2AE 



1 + 



A 



{A 



Csl) 



(4) 



where (Zcsi) = 54 is the average atomic number of the Csl medium. A^Ruth, 
the simple, analytical alternative to the volume concentration oc Sn, is 
plotted in Fig. 2b) for Wd = leV; the different atomic number curves are 
practically superposed - dashed hne. In the same figure, the first term of the 
above equation oc AZ^/E, weighted also by 1/Z^, was represented with a 
dotted line. This result provides a good approximation to A^Ruth- 



3.4 Calculation of 



The discrete functions in Fig. 1 can not be explained in terms of the carrier 
diffusion process only. Let us suppose that there is a carrier concentration 
threshold, below which quenching does not exist. The solution of a diffusion 
equation in cylindrical geometry predicts a concentration which diminishes 
with time and with the radial distance, but remains proportional to the initial 
density of carriers per unit length of the ion track. This fact may be only 
compatible with a monotonic decrease of the light output with the stopping 
power and not with the discrete functions mentioned above. In addition, the 
transverse outfiow of the carriers is somehow inhibited by space charge effects. 
A very high carrier concentration is maintained in the proximity of the incident 
particle path [47]. Consequently, the diffusion is not taken into consideration 
in the next calculation. On the other hand, we have seen in subsection 2.2 that 
this dependence is compatible with the knock-on electron generation. These 
electrons have ranges greater than r^, viz. energies higher than a cut off value 
Tcut, taking with them a fraction J-{E) of the energy deposited in the path 
element dx. 

If f3 and 7 are the kinematic variables of the incident particle, c is the speed 
of light in vacuum and u the unified atomic mass unit, the maximum energy 
Tmax which may be transferred to a stationary unbound electron of mass rUe 
may be calculated as T^ax — 2meC^/3^7^ in the "low energy" approximation 
2^me/{A X u) <^ 1 [48]. The condition that such an electron escapes the 
primary column and can be considered as a 5 - ray is that this energy is at 
least equal to Tcut', otherwise said: 

Tmax > = 2mecV|7|, (5) 
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Fig. 2. Specific energy loss divided by for different ions {AZ = 10) in Csl(Tl), 
versus the energy per nucleon. a) Specific electronic stopping power as in Refs. 
[44,45] - solid lines; approximate specific electronic stopping power used in calcula- 
tion oi - dashed lines, b) Specific nuclear stopping power as in Ref. [44] - solid 
lines. The dashed curves concerning the density of lattice ions displaced by Coulomb 
scattering (Eq. (4) for Wd = leV) are indistinguishable for all ions, and well ap- 
proximated by the first term of the mentioned formula - dotted line. For a) and b), 
see subsection 3.3 for explanations. 

where the index 5 at the kinematic variables refers to the minimum energy per 
nucleon eg of the incident particle that still generates 6 - rays ~ 2e^/ (uc^) 
in the non-relativistic approximation). (— dE/ dx)T<Tcut) defined in the previ- 
ous subsection, provides the specific electronic stopping power deposited inside 
the primary column. The energy deposited by the S - rays outside the primary 
column is (-dE/dx)^ = (-dE/dx)e - (- dE/ dx)T<Teuf I^om Eq. (2), the 
factor T becomes in this case: 



(-dE/dx)5 1 ln(^)-/32 + /3^ 



2 Tc 



cut 



(- dE/ dx)e 2 i ln(2iri«r^3,) - /32 - f ' 



(6) 



(A being the density effect in Bethe-Bloch formula [48]), which in the nonrel- 
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Fig. 3. a) The fractional energy loss carried outside the primary column by the 6 
- rays; it depends on the fragment velocity oc ^/E/A, but not on its identity (Eq. 
(7)). Different curves correspond to different hypothetical values of Tc„j, the cut 
off value in the electron energy above which they may be considered as 6 - rays 
escaping the highly ionized fiducial cylinder of radius Tc along the fragment wake, 
b) Tcut - solid line - and rc - dashed line - against the fragment energy per nucleon 
threshold necessary to generate S - rays. 

ativistic approximation <C 1 reads: 



Eq. 7 is valuable only for (3 > (3s, under this limit is zero. The only unknown 
quantity in the expression of ^ is eg, which will be one of the free parameters 
of the model. The factor vs E/A is plotted in Fig. 3a) for different threshold 
energy per nucleon es values. 




(7) 



2ln(^/3|)+ln(g) 
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3.5 Estimation of Tc 

One may calculate, from equation (5), the minimum 5 - ray energy values Tcut 
required to escape the primary column, as a function of e^. In the knock-on 
electron energy interval of interest here, it is found experimentally that the 
practical [12] range - energy relation can be reasonably described in various 
stopping media by a function of the form: i?e = hT"' [50] , where both b and n 
are constants [51,52]. With in mg cm~^ and T in keV, and considering Csl 
as having an intermediate mean atomic number between those of aluminium 
and gold, the two extreme media treated in Refs. [51,52], we find the particular 
value b = 0.016 by keeping n=1.35 [12,51] as universal power index for all the 
stopping media. Consequently, and without taking into account the angular 
distribution of the scattered 8 - rays, one may estimate the corresponding 
primary column radius r^- Both quantities: Tcut and Tc are plotted in Fig. 3b) 
versus e^. 

3. 6 The formalism 

Inside the primary column, the initial (t=0) volume concentrations of the 
created pair carriers: electrons (e) and holes (h), at the coordinate x along 
the particle path, are proportional to the electronic stopping power per unit 
length pSf.: 



for holes and electrons, respectively. The number of knock-on electrons per unit 
length escaping out of the primary column (as long as E /A > es) is TpS^/e (e 
being the mean ionization energy). They will produce other ionizations outside 
the primary column. The radial diffusion is not treated here. 

The basic assumptions concerning the thallium activator centres in the cae- 
sium iodide crystal is that they may exist in two different states: one acting 
as hole trap - of concentration AiV^ - and the other one - of concentration Na 
- which acts as electron trap, to produce afterwards the scintillation. 

The excited system inside the primary column has the tendency to come 
back to the equilibrium state. The electrons may be trapped by thallium 
electron traps with a probability AabA^a, or by defects with a probability 
Aoe (-^D + -^n), or they may directly recombine with holes at ionized cation 





ne(a;, 0) = neo(a;) = (1 - ^{x)) 



(1 - J^{x))no{x) 



(9) 
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and anion sites of the lattice, with a probabihty A^nh {x, t). The holes, in their 
turn, may be trapped by the thallium hole traps with a probability AakAA^a, 
or by defects, with a probability A^h (A^d + A^n) or they may recombine with 
the electrons, with a probability A^ne{x,t). The equations concerning the 
time variation of the carrier concentrations inside the primary column read: 



^^^^^ = [AAeA^A + Adc (A^d + Ar„)] n,{x, t) 

+ A^_nh{x,t)ne{x,t) (10) 



- i^^i^ = [AAhAA^A + Aoh (A^D + Ar„)] n^x, t) 

+ ARne nh(x,i). (11) 

By disregarding the activator centre depletion [24], the coefficients AacA^a and 
AalAA^a niay be considered as being approximately constant. For the sake of 
simplicity, we make the same hypothesis for Aoe (Ad + An) and Aoh (Ad + An)- 
In this way, we must solve two coupled first order differential equations with 
constant coefficients. The initial conditions are given by the expressions (8,9), 
where > for E/A > es and J-' = for E/A < e^. Of course the number 
of defects is much smaller than that of the activator centres, but we shall see 
in the next section that the nuclear induced defects play a role in describing 
the total light output at lower incident energies. In order to keep the number 
of free parameters in the resulting fit procedure as low as possible, we make 
the hypotheses that: Adh ~ Adc = Ad- Our tests have shown that the ratio 
(AAhAA^A + ADAD)/(AAeA^A + ADAD) ~ 1 and that the shape of the total light 
output is not critically sensitive to its magnitude around this value. Thus, in- 
stead of keeping this ratio as a free fit parameter, we fix its value at 1. We 
make the notations: Aacdr — Aac,d.r/(AacA^a + AdAd). Under the above 
assumptions, one may calculate the first derivative of the hole concentration 
with respect to the electron concentration: 

duh (x, t) ^ (1 + A'^N^ + A^rUc (x, t)) nh (x, t) 
dne(x,t) {l + A'^Nn + A'^ni,{x,t))n^{x,t)' 



The carrier concentrations are correlated. At every moment t and position 
X, the hole concentration rih {x, t) may be considered as a function of the 
electron concentration {x, t), expressed by a Taylor series expansion around 
the initial value at the x coordinate rie {x, 0) = rieo . By keeping only the zero 
and first order terms in the expansion, viz. by considering that ni^{x,t) is 
linearly dependent on 11^ {x, t) and by replacing it in the rate equation for 
the electron concentration Eq. (10), one may determine the time dependent 
solution: 
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ne{x,t)^ (13) 

p{x) {1 — (,t)) no (x) 
(1 + A'j^N^ (x) + A'^Tio (x)) exp (p (x) (Aas^Va + Ad^Vd) t) - y (x) A'rUq {x) 

with p(a;) = 1 + AJ-,A^n(a;) + (1 — y{x)) A'^no{x) in every point x of the ionizing 
particle trajectory, and y{x) = 1 — (x) A^no(x)/ (1 + A[)N-a{x) + A^no{x)). 

For the 6 - rays, which are supposed to produce new ionizations outside the pri- 
mary column, the new created carrier concentration equations are completely 
decoupled for electrons and holes, because there is no radiation damage nor 
are there recombination processes. In the electron equation of interest, only 
the first two linear terms (concerning the activator and permanent defect sites) 
on the right side of the Eq. (10) have to be kept. Thus, the solution will be a 
simple exponential. 



3.7 The light output 



The double differential hght output ( d^L(x, t)/ dxdt)p in the primary col- 
umn is given, within a multiplicative constant "a" (related to the energy - 
light unit conversion, the light collection in connection to the shape of the 
crystal and the PMT gain), by the first of the linear terms in equation (10): 
ae7rr'^A\eNAne{x,t). Integrated over time between and oo, this term leads 
to the differential light output in the infinitesimal element dx of the primary 
column: 

1 + A'^N^jx) + A'^nojx) 
1 + A'^N^{x) + (1 - y{x)) A'^noix) 



The differential light output generated by the 5 - rays in the same infinitesimal 
element dx (determined by integrating over time the corresponding exponen- 
tial solution of electron rate equation) is: 




= aA'j^^N\J^{x)pSe — ae-KrlA'j^^NAJ^{x)no{x). (15) 

s 



The total differential light output ( dL(x) / dx) is the sum of the light produced 
inside the primary column and outside it. The integration over the range of 
the particle R{Eq) or, by changing the variable, over the energy Eq leads to 
the total light output expression: 
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anSe{E) \ ^ dE 

l + a^S^iE) + anS,iE)) 1 + |1| 

y{E)anSe{E) dE 
1 + a,S^{E) + a^S,{E) J 1 + f||} 

(16) 



where y{E) = 1 - J^{E)aKSe{E)/{l + a^Sr^iE) + aKSe{E)). This hght output 
expression contains four model parameters: = aA'p^^N^p, cir = A'j^p/ (evrr^), 
On = Aqp/ (evrr^) x {e/e^), written Ukc this because Nn{E) = pS-a{E)/ (enVrr^) 
(here, en is the mean energy required to produce a dislocation) and {Eg — 
A X e^). These parameters denote meaningful quantities, og is the "gain" 
parameter, related to the energy-light unit conversion, light collection, con- 
version factor for producing photoelectrons and the electronic gain, or, are 
the "recombination quenching" and "nuclear quenching" parameters respec- 
tively, related to the processes which divert part of the deposited energy - 
what we call "quenching" . is the energy per nucleon above which the 5 - 
rays play a role in the scintillation mechanism. 

The first term in Eq. (16) refers to the last part of the particle range, where 
the energy has decreased below the b - ray production threshold E < Av. 
and thus T{E') — 0. The last two terms concern the first part of the trajectory 
(i? > A X and T{E^ > 0) as follows: the second term reports on the light 
output originating inside the primary column, while the third one is connected 
to the light produced by the knock-on electrons outside the primary column. 
For the low incident energies Eq < A x es, T{E) — along the entire range 
of the particle and only the first term will be present. 



4 Compcirison to the experimental data 

4-1 The light output - energy relation 

It is possible now to confront the light output - energy relation in Eq. (16) to 
experimental dependence, provided that the above four coefficients are known. 
Actually, they are determined, from data concerning ions with Z < 45, as fit 
parameters by a minimization procedure using the MINUIT package from 
CERN library. The experimental light output values Qo) and the correspond- 
ing energy Eq deposited in a CsI(Tl) crystal of INDRA (see the accompanying 
paper [1]) are used for the comparison. The integral in expression (16) is per- 



L = aG 



/; 



aKSe{E) 



In 1 



+ 



Es 

Eo 



is' 



1 - J^jE) 
y{E)anS,{E) 

T{E) dE 



In 1 



+ 



Sn(E) 
5e(E) 
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Table 1 

Fit parameters: ac, or, e^, gain, recombination quenching, 5 - ray production 
energy per nucleon threshold and nuclear quenching, respectively, for module 2 of 
ring 3. The errors on the parameters (one unit on the last digit) are only statistical. 
For exact calculation (Eq. (16)): a) cx 5*11; b) A^n ~ A^Ruth and Wd = 1 eV, were 
assumed, c) Pure recombination case. For the first order approximation calculation 
(Eq. (18)): d) oc Sn- See subsections 3.7, 4.5 for explanations. 





a) 


b) 


c) 


d) 


OG [a.u.] 


20.00 


19.94 


20.57 


18.74 


Or [eV~^10^^atoms~^cm^] 


2.82x10-2 


2.81x10-2 


3.04x10-2 


1.26x10-2 


es [MeV/u] 


1.96 


1.95 


2.08 


1.79 


On [eV~^10^^atoms~^cm^] 


3.9x10-1 


4x10-1 [a.u.] 




3.9x10-1 



formed numerically. The resulting fit parameter values are given in Table la). 
The fourth parameter an was held constant after a preliminary analysis. Us- 
ing these values in expression (16), the total light output was recalculated and 
compared to the experimental Qq values; the results are plotted versus Eq in 
Fig. 4 - solid lines. There is excellent agreement, within 3% for most data, as 
will be shown later. In the inset of figure, the low energy results are shown 
and compared to the case where the nuclear quenching term a^Sa{E) is disre- 
garded (dashed line) . The small unrealistic bump of this pure recombination 
quenching case, induced by the Bragg peak in the specific electronic stopping 
power Sc{E), is completely eliminated by including a^S^^{E). Thus, the nuclear 
quenching term may be considered as a "fine tuning" device allowing a more 
accurate description of the experimental light output at the lowest energies. 

Finally, when fixing at the proposed value (3.9x lO^i eV^il0i^atoms^icm2), 
it is very important to note that light ions up to Z = 7 are necessary but 
sufficient to determine the gain parameter qq to within 10% and the and 
eg parameters to within 25% of the above mentioned values. If three of the 
parameters: a^, es and are fixed at the values listed in Table la), the 
gain parameter aQ is determined to within 2.5% when the experimental data 
employed in the fit procedure are restricted to light ions Z < 7 too. 

When the defect concentration created by nuclear collisions, instead of being 
taken as oc Sn{E), is calculated from the Rutherford diffusion probability - 
Eq. (4) - the fit parameter values are practically the same (eventually except 
for an because iVRuth may be used in arbitrary units). These values are shown 
in Table lb) for Wd =1 eV. The hght output values calculated in this case 
are practically identical to those from the previous recipe. Thus, we conclude 
that the simple expression of iVRutiu very well approximated by its first term 
oc AZ'^/E as shown in Fig. 2b), may be used for the nuclear quenching term 
estimation. In this case, the value of the fourth parameter a„ depends on 
the value used for Wd' the smaller Wd is, the higher the number of created 
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imperfections (defect or large amplitude vibration levels) and the smaller the 
determined fit parameter value, and vice-versa. These two quantities cannot 
be kept simultaneously as free parameters in the fit procedure and one has to 
choose, for example, a value for Wo- 

In the end, we have to note that the agreement of calculation with data is 
rather good (except for the lowest energies) even when the nuclear term is 
disregarded, with the advantage that the number of fit parameters reduces to 
three. The resulting parameter values are also presented in the Table Ic). 

The quality of the above approach, which will be termed the recombination 
and nuclear quenching model (RNQM) in the following discussion, is illus- 
trated also in the next subsections by three direct consequences. 

4.2 The 5 - ray effect 

The fraction of the deposited energy J^{E) transferred to 5 - rays increases 
with E/A - Fig. 3a), and has the same value irrespective of Z. For two differ- 
ent fragments having the same average stopping power oc AZ'^/Eq, the higher 
Z fragment has a higher energy per nucleon, a higher J-'{E) value and, conse- 
quently, a higher experimental efficiency Qq/Eq, as shown by symbols in Fig. 
1. The calculated scintillation efficiency L/Eq provided by RNQM shows the 
same behaviour. The total light output is given by expression (16), with the 
nuclear quenching term a^S^{E) or a,^NYtMt\i{E) (for W-q = 1 eV - 1 keV). 
The results are shown by the curves in the upper box of Fig. 1 for a few of 
the heavier fragments, where the effect is quite pronounced. As long as the 
average specific electronic stopping power is not too high, they qualitatively 
follow the experimental data. In the latter domain, the curves are increasing 
when the nuclear quenching term is neglected. 

4-3 Reaction product identification in a Si - CsI(Tl) map 

At forward angles (rings 2-9), the INDRA array is composed of three detec- 
tion layers: ionisation chamber, 300 /xm thick silicon detector (Si) and Csl(Tl) 
scintillator. The reaction products cover a large range in energy and atomic 
number and they can be well identified in a Aii^gi ~ Qo map - symbols in Fig. 
5 - where AE'si is the energy lost in the silicon detector. With the values of 
the fit parameters aG,aR,an,e5 in column a) of Table 1 and the expression (16) 
for the light output, a Ai?si — L map was calculated - solid lines in Fig. 5. 
The cases b), c) of Table 1 give results practically identical to the ). The 

calculated curves have excellent agreement with the experimental loci, within 
AZ f« 1 for Z 50. 
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Incident Particle Energy Eg (l\/leV) 

Fig. 4. Total light output Qq against the initial particle energy £'0 for different ions 
(AZ = 5): the symbols are experimental data from the system Xe + Sn at 32 and 50 
AMeV [1], the solid lines are exact calculation (Eq. (16)) and the dotted lines are first 
order approximation calculation (Eq. (18)). Inset: zoom at low energies; comparison 
to the case where the nuclear quenching term would be neglected (dashed line) for 
Z = 30. 

4-4 Energy determination 

Having determined the four fit parameters, we may recalculate by Newton 
iterative method starting from Qq, the incident energy -Ecaic of the reaction 
product, to be compared to the true energy Eq. For most of the data, the 
relative deviation of the calculated energy versus the true energy per nucleon 
is less than 3% - as shown in Fig. 6a) for several ions with 5 < Z < 50. Only 
at low energies E < 5 AMeV, the discrepancy reaches rather large values: 5% 
- 10%, and up to 15% in the worst cases - see Fig. 6b). Note that this roughly 
corresponds to the true energy accuracy, shown by the regions between the 
two solid lines in Fig. 6a) [1]. 

These results are arguments in favour of the present "quenching" model which 
accounts for the main experimental evidence concerning the scintillation yield 
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Light Output Qg (a.u.) 

Fig. 5. A AEsi — Qo map (module 2, ring 3 of INDRA) from the sytem Xe + Sn 
at 32 and 50 AMeV. Conventions of presentation are the same as in Fig. 4. 

of the CsI(Tl) crystals in the yellow band of interest. The ratio an/a^ could 
give an idea of the weights of the two processes: recombination and nuclear 
quenching, provided that the ratio of the mean energy required to displace 
a lattice ion from its site to the mean ionization energy en/e is known. The 
parameter ~ 2 MeV/u, found by the fit procedure, to which corresponds a 
minimum energy of 4.4 keV for the 6 - rays, allows to roughly estimate - as pro- 
posed in subsection 3.5 - the radius of the primary column ^260 nm. This 
value is one order of magnitude greater than it has been found in Ref. [12]. 
Ref. [47] estimates the same quantity in Si to be 1.15/xm, with the connected 
S - ray energy ^ 10 keV. For electrons having this energy, the range in Csl is 
of about 770 nm, of the same order of magnitude as the value found in the 
present work. The results are given in Table 2. By means of the estimated 
value Tc ^260 nm, one may determine the order of magnitude of the carrier 
concentration uq generated inside the strongly ionized cylinder. For energies 5 
< E < 80 AMeV, one finds 10^^ < uq < 10^^ carriers cm'^ when 2 < Z < 50, 
i.e. Uq < Nao- This fact a posteriori supports disregarding activator centre 
depletion. 
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Fig. 6. Deviations (in %) of the energy values determined with Eq. (16) from the 
true energy values for several products of the reactions Xe + Sn at 32 and 50 AMeV: 
a) deviations plotted against the product energy per nucleon: symbols; the regions 
between curves show the accuracy of the true energy per nucleon (see subsection 
2.2 of the accompanying paper [1] for details); b) deviation histograms. 



Table 2 

Estimation of the radius Tc of the transverse area of the primary column: a) Ref. 
[12], b) prese nt work c) Ref. [47], 

a) Ref. [12] b) Present work c) Ref. [47]* 

rc [nm] 40 260 770 

* extracted from data in Si 



4-5 The first order approximation of the light output expression 



The quenching terms in Eqs. (14) for the differential light output and (16) 
for the total light output were evaluated by means of the parameter values 
from Table la) and stopping powers in the domain of our experimental data: 
On'S'n < OR-S'e <S 1. In the first order approximation of the logarithm, the 
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differential ligfit output in tfie primary column becomes: 

— ~«G7- — ^T^^-^^T^S,{x). 17 
\ dx / 1 + ar,S^{x) + aKSei^x) 



In the low energy case Eq < A^e^ for which J^{x) = 0, the above expression is, 
except the nuclear term, identical with the formula of Birks [7] from subsection 
2.1. In the same approximation, the total light output reads: 



L = aG 



Eg 



X 



dE 



1 + a^S^{E) + anS,{E) 1 + S^{E)/S,{E) 



+ 



Eo 

I 1 + a 



1 - J^{E) 



dE 



X 



Es 
Eo 



+ 



+ a^S^(E) + aj,S,(E) 1 + S^(E)/S,(E) 
T{E) dE 



Es 



1 + S^{E)/S,{E) 



(18) 



The resulting fit parameter values are given in Table Id). The quality of the 
results (total light output, fragment identification) obtained in the first order 
approximation is very close to that achieved by means of the exact treatment, 
as shown in Figs. 4, 5 - dotted lines. Formula (18), as well as expression (4) 
in the mentioned approximate form, will allow the analytical integration of 
the total light output, useful for the applications presented in section 4 of the 
accompanying paper [1]. 



5 Conclusions 



The nonlinear response of CsI(Tl) was quantified in terms of the competing 
deexcitation transition probabilities of the excited fiducial volume around the 
path of a strongly ionizing particle. A simple mathematical formalism was 
developed to determine total light output based on a numerically integrable, 
four-parameter dependent expression. The formalism has allowed the inclusion 
of energetic knock-on electrons. Consequently, their effect in the scintillation 
efficiency was quantitatively taken into account, in a fully consistent way, for 
the first time. By a fit procedure, the values of the four parameters were deter- 
mined. Of physical interest is the energy per nucleon threshold for producing 
5 - rays by the incident particle, = 2 MeV/w. This result also allows the 
rough estimation of the primary column radius rc=260 nm along the wake 
of the ionizing particle. The predictive character of the model was verified by 
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using telescope-type {AEsi — Qo) reaction product loci data, the agreement 
calculation - experimental data being within AZ ~ 1 for Z ^ 50, as well as by 
using deposited energy data, in the range 1-80 AMeV. A comparison of the 
model predictions to the latter experimental data resulted in variations of the 
order of the CsI(Tl) scintillator resolution 3%, and within 10% - 15% when 
the reference energy was known with only a comparable, poorer precision. The 
effect of the S - rays in the scintillation efficiency as a function of the average 
stopping power in the entire range is quite well described too. In fact, the 
knock-on electrons play an important role in the light output, especially for 
heavy reaction products. In the accompanying paper [1], an approximate but 
analytically integrablc formula will be derived from the proposed formalism. 
Up to now, it was used for rapid calibration and identification procedures in 
the INDRA An array. In the future, in relation with computer capabilities, we 
will use the exact formula proposed here. 
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